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Abstract. At a conference devoted to ever deeper surveys hunting for ever 
more distant galaxies, I posed a question for which a concensus view has been 
difficult to reach: 'Is there evidence for, or can we rule out, a significant popu- 
lation of dust-obscured starbursts at z > 3?' If, as seems likely, submm-selected 
galaxies are proto-ellipticals, one of the biggest unanswered questions is whether 
a significant fraction form at very high redshift - perhaps by the collapse of 
single gas clouds - or whether the entire population forms over a range of red- 
shifts, primarily via merging. The latter scenario is favoured strongly by existing 
data for the majority of bright submm galaxies - mergers are common; typical 
spectroscopic redshifts are in the range 1-3. However, our reliance on radio 
imaging to pinpoint submm galaxies leaves open the possibility of a significant 
population of very distant, massive starbursts. To rule out such a scenario re- 
quires a completely unbiased redshift distribution for submm-selected galaxies 
and this is unlikely to be forthcoming using conventional optical/infrared spec- 
troscopic techniques. Here, I summarise recent attempts to close that door, or 
pass through to an early Universe inhabited by a significant population of col- 
losal, dust-obscured starbursts. I conclude that the door is only barely ajar; 
however, the idea of galaxies forming, near instantaneously, on a collosal scale is 
not dead: I show an SMG imaged by MERLIN+VLA for >lMs - the deepest 
high-resolution radio image so far obtained. The data reveal a galaxy-wide star- 
burst covering >10kpc. Thus interpreting SMGs in terms of compact ULIRG- 
like events may not always be appropriate, as one might expect when considering 
the Eddington limit for >1O 3 -M yr -1 starbursts. 



1. Introduction 

Surveys with bolometer arrays at mm and submm wavelengths are potentially 
sensitive to dusty objects at extreme redshifts, galaxies that drop out of sur- 
veys at shorter and longer wavelengths due to obscuration and unfavourable K 
corrections. The first cosmological surveys, using the SCUBA and MAMBO 
cameras, radically changed the accepted picture of galaxy formation and evolu- 
tion, moving away from the optocentric view of the last century. The discovery 
of so-called 'SCUBA galaxies' (Smail, Ivison &; Blain 1997), now better known 
as submm galaxies or SMGs, was greeted with surprise due to the remarkable 
evolution in dusty, starburst galaxies implied by such a large source density at 
the flux levels accessible to the first generation of bolometer arrays. 

Sadly, SCUBA has been retired. SCUBA-2 will revolutionise submm as- 
tronomy once again, but not until the winter of 2007-08. At the depths that 
were accessible to SCUBA, and to cameras like MAMBO-117, SMGs constitute 
around half of the submm background and thus at least a quarter of the total ex- 
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Figure 1. 850- fim images of distant, optically-selected quasars (Priddey, 
Ivison & Isaak, in preparation). Left: SDSS J075618.13+410408.5 at z = 5.1; 
right: SDSS J104433.04-012502.2 at z = 5.8, with a nearby SMG that may lie 
at the same redshift (see also Ivison et al. 2000; Stevens et al. 2003). 



tragalactic background. The relation between metal production and background 
radiation then implies that roughly 25% of all the stars ever formed have formed 
in extreme objects like this (Eales et al. 1999). The combination of this argu- 
ment, their space density and their high bolometric luminosities made it hard 
to resist the conclusion that these objects are ellipticals seen during a phase in 
which a significant fraction of their stars are forming, loosely 'proto-ellipticals'. 

Initial excitement was tempered by the first efforts to study SMGs at optical 
and infrared wavelengths. Early reports, backed up with a study in the GOODS 
North field by Hughes et al. (1998), suggested that the majority of the submm 
population had no plausible counterparts at magnitudes routinely accessible to 
optical surveys. Attention was diverted to various redshift engines and broad- 
band photometric techniques (e.g. Townsend et al. 2001; Wiklind 2003). As a 
result, only a handful of detailed studies were attempted, often for extreme and 
possibly unrepresentative galaxies. 

Progress resumed when deep, wide-field, high-resolution radio imaging of 
submm survey fields were shown to recover roughly half of the SMGs in each 
field, with an astrometric precision of ~0.3" (e.g. Ivison et al. 1998, 2002). 
Combined with the submm flux density, radio data provide a rough estimate of 
redshift (Carilli &z Yun 1999). They also enable refinement of submm samples, 
increasing the detection fraction to two thirds of SMGs at 850-/im flux density 
levels in excess of ~5mJy (Ivison et al. 2002; Pope et al. 2006). Accurate posi- 
tions in hand, bright SMGs were found to be a diverse population — some are 
quasar- like, with broad lines and X-ray emission (e.g. Knudsen et al. 2003), some 
resemble BAL quasars (e.g. Vernet & Cimatti 2001), some are morphologically 
complex (e.g. Smail, Smith & Ivison 2005), some are extremely red (e.g. Smail 
et al. 1999; Webb et al. 2003; Dunlop et al. 2004); some bear the signatures of 
obscured active nuclei and/or powerful winds (e.g. Smail et al. 2003). 
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Figure 2. Possible spectral energy distributions (SEDs) describing the emis- 
sion from a typical SMG at z = 2.4, with 850-/im and 1.4-GHz flux densities 
of 6mJy and 50/iJy, respectively, taken from C05. Overlaid are SEDs for 
Tdust = 25, 32 and 45 K. Upper panel: SEDs normalised to the radio datum, 
emphasising how galaxies with hotter dust (and lower implied dust masses) 
are missed in the submm at z > 2. Lower panel: SEDs normalised at 850 fim, 
highlighting how sources with cooler dust are undetectable in the radio at 
high redshifts. The need to pinpoint the counterpart of an SMG means that 
the coldest SMGs, which will lie below typical radio flux limits, will be missed 
by spectroscopic surveys. This bias becomes potentially serious at the highest 
redshifts where even relatively warm SMGs fall beneath typical radio limits. 

Not surprisingly, given their median magnitude (it! = 25, Vega), spectro- 
scopic redshifts have been difficult to determine. The most ambitious survey 
of SMGs was undertaken by Chapman et al. (2003, 2005 — hereafter C03, 
C05). The median redshift was ~2.2 for S'sso/im > 5-mJy galaxies selected us- 
ing SCUBA and pinpointed at 1.4 GHz or, in some cases, selected as optically 
faint radio galaxies and detected subsequently with SCUBA. The accurate red- 
shifts reported by C03 and C05 facilitated the first systematic measurements of 
molecular gas mass for SMGs (10 10 -10 11 M ) via observations of CO (Neri et 
al. 2003; Greve et al. 2005; see also early work by Frayer et al. 1998, 1999), as 
well as constraints on gas reservoir size and dynamical mass (Tacconi et al. 2006 
- though see later). The data suggest SMGs are massive systems and provide 
some of the strongest tests of galaxy-formation models to date. 



2. Current state of play 

In spite of this progress, a detailed understanding of SMGs remains a distant 
goal that many believe is now best achieved by detailed study of individual 
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Figure 3. The 850-/j,m/1.4-GHz flux density ratio versus redshift for the 
radio-detected SMGs with spectroscopic redshifts, taken from C05. The pre- 
dicted variation in flux density ratio for three SEDs are shown, including that 
of Arp 220. If the SMG population displayed a narrow range of SEDs we would 
see them tracing just such a track in this plot. The trend of submm/radio 
flux density ratio with redshift is flatter than the template tracks. The wide 
scatter, even at fixed redshift, suggests a wide variety of SEDs. The shaded 
region shows the ±<r envelope of the dispersion calculated in three redshift 
bins containing equal numbers of SMGs. The radio selection criterion for ex- 
isting spectroscopic work biases those surveys to galaxies with enhanced radio 
emission, such as Mrk231. A purely submm-sclcctcd sample would show an 
even wider range in submm/radio flux density ratio than the one shown here. 
The range of SMG SEDs means that simple photometric redshift techniques 
struggle to predict accurate redshifts for individual galaxies. 



objects. Confusion currently limits our investigations to the brightest SMGs 
(although surveys through lensing clusters have provided a handful of sources 
typical of the faint population that dominates the cosmic background, e.g. Smail 
et al. 2002; Kneib et al. 2004; Borys et al. 2004). We must also recall that 
selection biases have potentially skewed our understanding. Around half of all 
known SMGs remain undetected in the radio due partly, at least, to the lack of 
sufficiently deep radio data, which do not benefit from the same K correction as 
submm data. Although a handful of radio-undetected SMGs have been targeted 
spectroscopically, few have yielded reliable redshifts. Finally, there has been only 
limited coverage at red and IR wavelengths in existing spectroscopic surveys. 

Here, at a conference devoted to ever deeper surveys hunting for ever more 
distant galaxies, I posed the question: "Is there evidence for, or can we rule out, 
a significant population of dust-obscured starbursts at z > 3?" 

It is difficult to present a concensus view: this issue divides the submm 
community, with opinions correlated distinctly with age. Deep imaging was 
not possible in the pre-SCUBA era so submm cosmologists spent their time 
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hunting for thermal emission from distant quasars and radio galaxies (Fig. 1; 
e.g. McMahon et al. 1994; Ivison 1995; Hughes, Rawlings & Dunlop 1997). In 
the eyes of the SCUBA generation, these saggy, old astronomers have found it 
difficult to accept that their dusty high-redshift AGN were rare freaks which 
modellers can label 'unrepresentative'. The C05 spectroscopic survey suggests 
that the vast majority of SMGs lie at z = 1 — 3 with a mere smattering at z < 1 
and z > 4: 20% at most (cf. Dunlop 2001). But is this necessarily so? Would 
we be aware of a significant population of distant, dusty starbursts, maybe at 
the >25% level, hiding within current SMG samples? 

Several potentially important effects - not least our reliance on radio imag- 
ing to pinpoint SMGs (see Fig. 2) - leave the door open for a significant pop- 
ulation of distant, massive starbursts. To rule out such a scenario requires an 
unbiased redshift distribution for SMGs and this is unlikely to be forthcom- 
ing by conventional spectroscopic techniques (though see Pope et al. 2006 for 
the best available attempt using conventional optical/IR photometric methods). 
Before ALMA, JWST or the ELT, and without the construction of dedicated 
wide-bandwidth cm/mm/submm spectrometers designed to solve this specific 
problem, the only way to make progress is to use the distinctive SEDs of high- 
redshift dusty galaxies as a method of estimating redshifts. 

3. Recent work 

Carilli & Yun (1999) and Eales et al. (2003) were amongst the first to address 
this issue, along with Hughes et al. (2002), Wiklind (2003) and Aretxaga et al. 
(2003), though cf. Blain, Barnard & Chapman (2003) and C05. The approach 
generally employs long- wavelength flux ratios as redshift indicators. Carilli & 
Yun pioneered the use of the S^o^m/Si^GHz flux density ratio (Fig. 3), men- 
tioned first by Hughes et al. (1998) and used later by Ivison et al. (2002) to 
determine a median redshift for SMGs of >2; however, this particular ratio 
flattens at higher redshifts and the flux ratio is not very useful at z > 2. 

Eales et al. (2003) developed a different redshift indicator: the ratio of 850- 
to 1200-/im flux density, as traced by MAMBO and SCUBA, respectively. This 
has complementary properties to the radio/submm flux ratio: it is useless as a 
redshift indicator at z < 2.5, but it should very good at higher redshifts. For 
answering the specific question of whether there is a significant population of 
dusty, starburst galaxies at z > 3, as one would expect if part of the elliptical 
population forms by monolithic collapse, this is the most useful redshift estima- 
tor. Fig. 4 shows the predicted 850-/1200- fim flux density ratio as a function of 
redshift for several template SEDs. Errors for individual redshift estimates are 
clearly very large, to put it mildly, but the error on the mean redshift of a large 
sample is much smaller. 

Eales et al. carried out SCUBA 850-/um photometry-mode observations (in 
which a single bolometer in each array is pointed at the target) of the brightest 
sources from one of the earliest blank-field 1200-^m MAMBO surveys. Using 
SCUBA to follow up MAMBO sources exploits two advantages over MAMBO 
follow-up of SCUBA sources: first, a 1200-//m survey is likely to contain a larger 
fraction of very-high-redshift sources than an 850-/um survey; second, the IRAM 
30-m MRT's resolution at 1200 ^m is 30% better than the SCUBA resolution 
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Figure 4. Predicted 850-/1200-/im flux ratio as a function of redshift from 
Eales et al. (2003). Asterisks show measured values of the flux ratio for 
MAMBO-selected sources, plotted arbitrarily at z = 0.2; dots: high-redshift 
quasars; horizontal line: range of spectroscopic redshifts for SMGs, circa 2003; 
thin lines: flux ratio versus redshift predicted for single-temperature dust 
models; dotted lines: (3 = 1; dot-dashed lines: (3 = 2. For both sets, the 
lowest and highest lines are for Td us t = 20 and 70 K, respectively, with other 
lines at intervals of 10 K. Thick lines: predictions based on the models of 
Dunne & Eales (2001). Continuous line: median predicted value of the flux 
ratio; dashed lines: lowest and highest predicted value at each redshift. 

at 850 fim which is crucial for accurate photometry. More accurate positions 
from radio imaging would be better still, but this would re-introduce our most 
pernicious bias. 

Eales et al. concluded that the 850-/1200-^m flux density ratios for 15 of 
their sample of 23 SMGs were much lower than expected for low-redshift galaxies 
and that these SMGs must either be at higher redshifts, z > 3 in fact, or have 
different rest-frame SEDs that local starbursts. Various systemic effects were 
investigated and ruled out, including concerns about absolute flux calibration, 
astrometric errors and flux boosting, the latter being a particular example of 
the process described by Eddington (1940) in which the statistical properties of 
a a distribution of experimental measurements are distorted by noise. 

Greve et al. (2004) took the next logical step, using unbiased blank-field 
surveys at 850 and 1200 /im (Fig. 5) to investigate the reliability of SMG samples, 
to analyse SMGs using flux ratios sensitive to redshift at z > 3, and to search 
for 'SCUBA drop-outs', i.e. z » 3 galaxies, where the far-IR bump has moved 
through SCUBA's 850-^m filter. For a sub-sample of 13 SMGs detected by 
both MAMBO and SCUBA, Greve et al. concluded that the distribution of 
850-/1200- j^m flux density ratios was consistent with the spectroscopic redshift 
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Figure 5. Left: MAMBO 1200-^im image of the Lockman Hole East field 
from Greve et al. (2004). Sources detected at >3.5<r are circled. The large 
white circle indicates the brightest MAMBO source that lacks a counterpart 
in the SCUBA 850-^m image of the same region: LH 1200.02, the second most 
significant source in the MAMBO image - see also Fig. 6. The 850-/im image 
is shown, at roughly the same scale, to the right, this time smoothed with a 
beam-size Gaussian (14.5" fwhm), from Scott et al. (2002). 




Figure 6. Spitzer imaging at 3.6, 4.5, 5.8, 8 and 24 /tm of the 'SCUBA 
drop-out', LH 1200.02 (courtesy: Eiichi Egami). The final panel shows a deep 
i?-band image. The counterpart to the SMG is far from obvious. 



distribution of radio-detected SMGs (C03, C05). For the 18 MAMBO sources 
not detected by SCUBA, the distribution of 850-/1200-/im flux density ratios 
was indistinguishable statistically from that of the SMGs identified by both 
MAMBO and SCUBA. 

Although Greve et al. (2004) disagreed fundamentally with the findings 
of Eales et al. (2003), it is easy to see that the sub-sample of SMGs detected 
by both MAMBO and SCUBA - which makes up less than half of the parent 
MAMBO sample - will (by definition) contain no 'SCUBA drop-outs'. At this 
stage I therefore consider that the door to an early Universe inhabited by a sig- 
nificant population of collosal, dust-obscured starbursts is still ajar. In fact, the 
sample contains at least one object detected robustly (6.8a) at 1200 fim and yet 
apparently absent in the SCUBA image - an object with the very characteristic 
that Greve et al. set out to find. This SMG (LH 1200.02 at 10 52 38.8, +57 23 
22, J2000) is circled in Fig. 5. It lacks both a robust radio detection and a con- 
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Figure 7. Postage stamps of galaxies from the dual-survey sample of Ivison 
et al. (2005). Extremely deep near-IR data are shown as a greyscale (K = 
20.9, 3cr, Vega), upon which 1.4-GHz contours are plotted. Open crosses: 
SCUBA positions; solid crosses: MAMBO positions; 8-arcsec-radius circles: 
average positions; boxes mark robust radio sources. A remarkable number 
of SMGs have counterparts that are barely visible at K ~ 21 (Vega). The 
median optical magnitude of this SMG sample is R — 25.0 (Vega). 



vincing Spitzer identification (Fig. 6). Elucidating its nature - hyper luminous, 
dusty starburst at z > 3 or an less distant object with significantly lower charac- 
teristic dust temperature, possibly <25K - has progressed at an excruciatingly 
slow pace. 

Ivison et al. (2005) presented a robust sample of bright SMGs selected using 
a dual-survey extraction technique, aiming to leave the regime where the modest 
significance of sources compromises the analyses via the presence of spurious 
sources or the effect of flux boosting. We would expect the resulting sample to 
be detected in well-matched radio imaging (o"850/im/o"i.4GHz ~ 500). The goal 
was to determine the true fraction of radio drop-outs amongst SMGs, as well as 
some practical information such as the intrinsic positional uncertainty for SMGs 
in the absence of radio/IR counterparts. 

The dual-survey selection technique yielded 19 SMGs of which a high frac- 
tion (~80%) had radio counterparts (Fig. 7). Given the likely ~10% contami- 
nation by spurious sources, this suggests that very distant SMGs (z S> 3) are 
unlikely to make up more than ~10% of the bright SMG population and that 
almost all of the Si mm > 4 mJy SMG population is amenable to study via the 
deepest current radio imaging, e.g. the <Si.4gh z = 2.7 /iJy r.m.s. imaging of the 
Lockman Hole field at 1046+59 by Owen et al. (in preparation). A number of 
caveats remain, however: for one, an AGN contribution to the radio flux density 
could mask a high-redshift population. 
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Figure 8. 5" x 5" MERLIN+VLA 1.4-GHz image (0.2" fwhm), with un- 
precedented sensitivity (>lMs of integration), of an SMG at z = 1.2 (Biggs 
& Ivison, in preparation). The radio emission, which should be thought of 
as a high-resolution surrogate for the far-IR emission, is extended (unlike the 
optical) on a scale of >2" (>15kpc), with several compact knots. 



4. SMGs as scaled-up ULIRGs? 

It is often argued that SMGs resemble scaled-up ULIRGs, with compact, intense 
starbursts occurring on scales of a few kpc or less. Measurements of small source 
sizes for their radio or CO emission are often cited to support this notion (e.g. 
Tacconi et al. 2006). Evidence exists, however, that this compact mode of galaxy 
formation is not necessarily ubiquitous amongst high-redshift galaxies, nor would 
we expect it to be when we consider the Eddington limit for >10 3 -MQyr -1 
starbursts. Ivison et al. (2000) and Stevens et al. (2003, 2004) presented images 
of distant, dusty AGN - radio galaxies and X-ray-absorbed quasars - showing 
rest-frame far-IR emission on scales of 10s of kpc, though whether that emission 
is due to smooth, galaxy-wide events or due to a number of intense starbursting 
clumps has so far been impossible to ascertain. 

Bright SMGs for which high-resolution submm/mm imaging have been 
sought to address this question, thus far in vain, include HzRG850.1 (Ivison 
et al. 2000). Now, however, exceptionally deep, high-resolution radio imaging 
with the MERLIN array (Biggs & Ivison, in preparation) - an unprecendented 
1-Ms integration, combined in the uv plane with 270 ks of VLA data using a 
new technique - means that we can investigate the scale of the radio (and thus 
far-IR) emission in normal SMGs. Fig. 8 shows one example: LH 1200.08 at 
z = 1.2. Its radio emission is extended on a scale of at least 2" (>15kpc), 
with several compact knots. As such, it bears only a passing resemblence to 
low-redshift ULIRGs. 

It is interesting to speculate about how this object would appear if it were 
at observed with the resolution and sensitivity more usually available to us and 
it lay at the median spectroscopic redshift of SMGs? The answer: it would 
resemble many of the existing radio detections - a faint, barely-resolved radio 
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source displaying very little evidence of the galaxy-wide starburst apparent in 
this deep MERLIN image. 
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